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^ Abstract. Convex polytopes are convex hulls of point sets in the n-dimensional 

O space E" that generalize 2-dimensional convex polygons and 3-dimensional convex 

\^ polyhedra. We concentrate on the class of rt-dimcnsional polytopes in E" called 

v^ sign permutation polytopes. We characterize sign permutation polytopes before 

^SJ relating their construction to constructions over the space of quantum density matrices. 

Finally, we consider the problem of state identification and show how sign permutation 

(~| polytopes may be useful in addressing issues of robustness. 

J9 PACS numbers: 02.40.Dr, 02.40.Ft 

'~~' 1. Introduction 

(N 

^ From the Platonic solids to n-cubes, n-cross-polytopes and regular n-simplices, convex 

■^ polytopes as geometric objects have been a source of study since antiquity. Indeed, 'the 

^ Universe cannot be read until we have learned the language in which it is written. 

[^ It is written in mathematical language, and the letters are triangles, circles and 

O other geometrical figures' - Galileo. A celebrated feature of these convex geometrical 

,__! figures has been the symmetry they possess. In micro-biology where objects use the 

j> efficiency of symmetry to propagate themselves to the study of genetics where symmetry 

k/H communicates genetic information, symmetry has been a fundamental concept of science 

H and has provided an understanding for many of Nature's phenomena. But perhaps the 

most important class of symmetric polytopes is the class of regular convex polytopes in 

Euclidean spaces. 

In this paper, we concern ourselves with the class of n-dimensional convex polytopes 
in E" called sign permutation polytopes. We show that regular convex polytopes 
in Euclidean space can be derived from such polytopes, and we show that their 
characterization can be admitted by the partial order of weak majorization. The outline 
of this paper is as follows. Section |2] provides an introduction to material from discrete 
geometry which serves as a basis for our study. Section |3] introduces a theorem of 
Rado (Rado (1952)) and motivates the issue of n-dimensional polytopes of degree n 
in the n-dimensional space E" of section |4j Section |5] relates the construction of sign 
permutation polytopes to the space of density matrices in quantum mechanics, before 
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finally illustrating how sign permutation polytopes may be useful in addressing the 
problem of robustness of state identification. 

2. Preliminaries 

2.1. Convex sets & combinations 

Let n G N and denote by E" the n-dimensional Euclidean space, with origin a, scalar 
product (■ , ■ ) and induced norm || ■ || . Let a = (ai, a2, ■ ■ ■ , an) G E" be an n-tuple of real 
numbers. We say a G E" is a linear combination of Oi, . . . , a^ G E" if a = YlT^i "^i^i f^^ 
suitable scalars Aj G M. If such scalars Aj G M, i = 1, . . . , m, exist with Xll^i -^j = 1 then 
we say a is an affine combination of ai, . . . , a^ G E". The set of points oi, . . . , a^ G E" 
is said to be afiinely independent if a linear combination Xll^li -^j'^« with Xli^i ^i = 
can only have the value a when Aj = for i = 1, . . . ,m. For A C E", the set of 
all affine combinations of points of A is called the affine hull of A and is denoted by 
aSA. If a = Yl^i -^i^i f*^^ non-negative scalars Aj G M and Yl^i ^i = 1 then we say 
a G E" is a convex combination of ai, . . . , a™ G E". For A C E", the set of all convex 
combinations of any finitely many points of A is called the convex hull of A and is 
denoted by convA. By i/a,/3, we mean the hyperplane given by the linear equation of 
the form (a, x) = (3. Consequently, the hyperplane Ha^js can be expressed as the set 
{x G E"|(a,x) = /?} where a G E^yjo} and /3 G M. A closed half-space of E" is given 
by the set {x G E"|(a, a;) > /?}, a G E"\{o} and /3 G M, where its boundary is the 
hyperplane {x G E"|(a, x) = /?}, a G E"\{o} and /3 G M. For convex sets A,_B C E" 
satisfying Afl-B = 0, the Hahn-Banach separation theorem (Bishop and Bridges (1985)) 
states that there exists a unit vector -u G E" and a scalar /3 G M such that for all a G A 
we have {u, a) > (3, while for all a G -B we have {u, a) < (3. 

2.2. The symmetric group & majorization 

Let Sn denote the symmetric group of degree n defined as the group of all permutations 
of the integers l,2,...,n. Denote by n{a) = (Q;,r(i) , • • • j ^^^(n)), the permutation of 
a G E" by TT G Sn- Let a* = (aJ, . . . , a*) G E" be the point whose components are those 
of a arranged in non-increasing order of magnitude; a* > . . . > «* and al = a7r(fc) , for 
k = 1, . . . , n and n G Sn- Let b* = (/J^ , . . . , /?*) G E" be defined analogously. If the 
relations 

k<l k<l 

hold, for 1 < I < n, with equality for / = n, we say a is majorized by b and write a -< b 
(Marshall and Olkin (1979)). For a,b E E", we have it that 

a + b^a* + b* (2) 

since there exists a permutation n E Sn such that an{i) + /37r(i) > 0:71(2) + f3n{2) > • • • > 
a^in) + I3n(n) and Y.k<ii^Hk) + (^^(k)) < Y.k<ii^k + Pk)^ for 1 < / < n with equality 
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for / = n. Consider the set of points ai,a2, ■ ■ ■ ,am G E"; a^ = (aji, ai2, • • • ,am), 
i = 1,2, . . . ,m. Following from ^, we have it that ai + a2 + - ■ ■ + am -< al + al + . . . + a'^. 
Finally, if Oj -< b, i = l,...,m, then X]i<m •^«'^« "^ ^ ^^^ non-negative scalars Aj and 

}-^i<m ^i ~ -'-• 

^. 3. Convex polytopes 

Convex polytopes are non-empty compact convex subsets of Euclidean space described 
by a finite point set. Although the concept of convexity is elementary, polytopes 
possessing this property yield geometric bodies with rich structure. Two-dimensional 
convex polytopes are called convex polygons. The equilateral triangle constructed by 
Euclid in his first proposition is such an example. Three-dimensional convex polytopes 
are called convex polyhedrons and examples include the Platonic solids (Coxeter (1969)). 
Interestingly, while the polyhedron has eluded definition over the centuries, the theory 
of convex polyhedrons and polytopes is well understood and contributes to practically 
significant areas of combinatorial optimization and computational geometry. To see 
this, consider the permutahedron defined by taking the convex hull of all vectors that 
are obtained by permuting the coordinates of the vector (1,2, ...,n). A feature of 
the permutahedron is that the vertices are identified with the permutations of the 
symmetric group of degree n in such a way that two vertices are connected if and 
only if the corresponding permutations differ by an adjacent transposition (Gaiha and 
Gupta (1977)). 

A convex polytope is the convex hull of a non-empty finite set A = 
{ai,a2, ■ ■ . ,am} C E". In particular, the V-polytope is the set of points describing 
such a polytope in terms of its vertices and is given as 

' m m \ 

J2 ^i^'i I Ai > 0, ^ Ai = 1 I . (3) 

=1 j=i J 

On the other hand, a convex polytope may be described as the bounded solutions set 
of a finite system of half-spaces in E". In this instance, we say such a polytope is a 
Ti-polytope and is described as 

{a; G E" I (a^, x) < A for 1 < i < m}, (4) 

for Cj G E" and /3j G M. A basic result of convex polytopes maintains that we may regard 
a convex polytope as a bounded solution set of finitely many closed half-spaces in E" or 
as the convex hull of a non-empty finite set A C E". Remarkably, while the descriptions 
of a convex polytope as that of a H-polytope or as a V-polytope are equivalent, the 
computational complexity associated with describing each differs in the extreme. For 
instance, the n-dimensional cross-polytope as a V-polytope is defined as the convex hull 
of the set of 2n points 

7„ = conv{7r(±l,0,...,0),7r G ^n}, (5) 

whereas, as a "H-polytope, the n-cross-polytope is described by 2"' half-spaces 

7„ = {x G E" I {ax) < 1}, (6) 
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where a runs through all vectors in {—1,1}". The surface of octahedron, a 3-cross- 
polytope, consists of six O-dimensional faces called vertices, twelve l-dimensional faces 
called edges and eight 2-diniensional faces called facets. Generally, the faces of an 
n-dimensional convex polytope have dimensions — 1,0, l,...,n with —1 denoting the 
dimension of the empty set with the face of dimension j being called a j-face. Finally, 
for A C E", a convex polytope is called a /c-polytope if dim conv A = k. This implies 
that some (A; + l)-subset oi A is afiinely independent but no such (A; + 2)-subset is afiinely 
independent (Ziegler (1995)). 

3. A note on Rado's theorem 

Rado (1952) characterized the convex hull of the set of all permutations of any real 
n-tuple in terms of the Hardy-Littlewood-Polya partial order relation (see Hardy et al 
(1934), p. 45) for real n-tuples. We now state Rado's theorem and discuss the nature of 
the convex set established by the theorem. 

Theorem 3.1 (Rado (1952)) Let a G E" and Va be the convex hull of the set of 
permutations of a. Then x &Va if and only if x -< a. 

The polytope of Rado (1952) is the {n — l)-dimensional permutahedron of degree n 
embedded in the n-dimensional space E". Since the vertices are obtained by permuting 
the integers (1,2, . . . ,n), this permutahedron of degree n lies entirely in an (n — 1)- 
dimensional hyperplane consisting of all points whose coordinates sum to the integer 
1 + 2 + ... + n = n{n + l)/2. For the application in quantum state identification, we 
concern ourselves with the construction of n-dimensional polytopes of degree n in the 
n-dimensional space E" that exploit the key feature of Rado's theorem. We now show 
that a necessary and sufficient condition for such polytopes requires a set of (n — 1)- 
dimensional bounding facets to contain points whose coordinate sums differ. 

Proposition 3.2 Let a,b & E". Let Va and Vb denote the convex hull of the set of 
permutations of a and b, respectively. Then Va ri Vb = ^ if and only if X]fc<n ^k 7^ 



To establish proposition |3.2[ we make use of the following results. Let A = 
convjoi, . . . , a^} C E" be a non-empty convex set. By the relative interior of A, relintA, 
we mean those points a G E" for which a = Y^iLi ^i^i "with Aj > 0, for i = 1, . . . ,k, and 

Er=i A. = 1. 

Lemma 3.3 (Griinbaum (2003)) Let A, B G E" be nonempty convex sets. Then A and 
B can be properly separated if and only if 

relintA n relintfi = 0. (7) 

Lemma 3.4 If a,b G E" such that J2k<n^k = J2k<nt^k, then Va and Vb can not be 
properly separated by a hyperplane Hua, 

H^^^ = {zeE^ I {z,u)=a}, (8) 
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Figure 1. Truncated octahedron: sign permutation polytopc with basis vector 
a = (0,1,2) centered at the origin. The truncated octahedron is an Archimedean 
solid consisting of 8 regular hexagonal facets, 6 square facets, 36 edges and 24 vertices. 



for M G E" \ {o}, and a G 



Proof. See Appendix A. Proof of proposition 3.2 



Let us suppose that X]fc<n'^'s ~ ^k<nl^k- By lemma 3.4, the convex hulls Va and 



Vb can not be properly separated. Consequently, by lemma 3^, we have Vafl Vj, 7^ and 
the implication follows. On the other hand, suppose Va r\Vb y^ 0- Hence there exists 
X G Va n Vb, and therefore numbers Aj > 0, /ij > satisfying J2i<m ^i = J2i<m' f^i = ^ 



such that X 



Y.i<m ^i(^i and X 



Y.i<ra' P^i^i- FoT I < U and X = (Xi, X2, . . . , x„) G E'^ 



we have it that 



k<l 



2J 2J -^«"«fc 
k<l i<m 



J] A, J] a,fe < J] A, J] a* = J] 



a 



ki 



(9) 



i<m 



k<l 



i<m 



k<l 



k<l 



with equality for / = n. Similarly, for / < n, we have it that 



k<l 



k<l i<m' 



i<m' 



k<l 



i<m' 



k<l 



k<l 



with equality for / = n. The implication follows, and this completes the proof. D 

While proposition 3.2 illustrates a necessary and sufficient criterion to construct 
an n-dimensional polytope if given an initial (n — l)-dimensional permutahedron, it 
suggests a generalization of Rado's theorem is not readily available since the inclusion 
of an arbitrary point does not preserve the majorization criterion. For example, consider 
a, 6 G E*^ and further consider the polytope described by taking the convex hull of Va 
and Vb- Let us denote such a polytope by convjVa, V;,}. For some x E Va and y G V^, 
the question of whether a point 2 G E" is contained in the polytope convjVa, V^} can 
be verified by showing that there exists a suitable AgM, 0<A<1, such that 
z = Xx + {1 — X)y. Consequently, we can demonstrate that z -< Aa + (1 — X)b, and, 
hence, z G convjVa, V;,} but only at the insistence of a suitably chosen scalar A G M. 
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4. Sign permutation polytopes 

Let us instead consider tlie task of constructing a polytope wliose vertices are not 
arbitrary ra-tuples in E". We adapt the theorem of Rado (1952) to construct an n- 
dimensional polytope of degree n in E"^ with vertices given by sign permutation points 
of E" and an interior that admits a characterization with respect to the partial order of 
weak majorization. 

Consider the convex polytope whose vertices consist of the set of sign permutations 
on the integers (1,2,..., n) defined as the set of all permutations of (±1, ±2, . . . , ±n). 
We call such a polytope consisting of 2"ra! vertices a sign permutation polytope, see 
figure [ij The sign permutation polytope is n-dimensional convex polytope of degree n 
in the n-dimensional space E". More generally, denote by V±^(a) the sign permutation 
polytope given as the convex hull of componentwise sign permutation changes of a G E"; 

V±7T{a) = conv{(±a^(i), ±a,r(2), • • • , ±a^(„)), ir G Sn}. (H) 

Let a* G E*^ denote again that vector whose components are those of a arranged in 
non-increasing order of magnitude. Let 6 G E" be defined analogously. If the relations 

E"^^E/5^ (12) 

k<l k<l 

hold, for 1 < I < n — 1, and 

E"^<E/5^' (13) 

k<n k<n 

we then say a is weakly majorized by b, and write a -<^ b. Let E" denote the set 
IE+ = {(ai,...,a„) I aj > 0, i = l,...,n}. 

Corollary 4.1 Let a,x E E" . Then x G V-i-7r(a) ^/ one? only if x -<^ a. 

Proof. If X G V±7r(a), then there are numbers Aj and vectors a^, i = 1,2, ... ,m, for 
m < 2"n! such that a, -K^ a, Aj > 0, J2i<m -^1 = 1 ^^^ ^ = X]j<-m '^i^i- Then since 
ai -<u) a, we have it that x = X]i<m ^i^^i "^^ ^i<m A^a = a for Ai > 0, z = 1, 2, . . . , m, 
and J2i<m^i ~ ^- ^^ ^^^ other hand, suppose x -<„, a while x ^ V±7r(a)- Using the 
Hahn-Banach separation argument of Rado (1952), the reverse implication is readily 
established, and the result follows. D 

We make note of the following relationship between the convex hulls of V±Tr(^a) and 

V±,r(z)- 

Corollary 4.2 Let a,x E E". Let V±Tr(^a) o,nd V±.,^i^x) denote the convex hull of sign 
permutations of a,x G E", respectively. Suppose x -<^ a. Then V±7r(z) ^ V±,r(a)- 

Proof. Given a,x E E" with x -<^ a, we have it that ±7r(a;) -<^„ x for tt G Sn- Noting 
that -<y^ is a partial order then transitivity of -<^ establishes ±7r(a;) -<w ±7r(a), and the 
result follows. D 

Remark 4.3 We note that Mirsky (1959) has constructed a polytope similar to the 
polytope of corollary \4 1\ restricted to the positive orthant of n-dimensional Euclidean 
space. Markus (1964) characterized the boundary of the sign permutation polytope. 
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A regular polytope generalizes the Platonic solids to arbitrary dimensions. For 
the n-dimensional spaces E" with n > 5, there are only three regular convex polytopes 
(Coxeter (1948)). These are the n-dimensional regular simplex, the n-dimensional cross- 
polytope and the n-dimensional cube. For a = (a, 0, . . . , 0) G E", the sign permutation 
polytope V±Tr(a) degenerates to the class of cross-polytopes /3„. The cross-polytope /3„ 
can be described in terms of the 2n vertex set {'7i{±a, 0, . . . , 0), tt G S'„}, or, equally, 
in terms of the 2" half spaces 7„ = {x G E" | (ax) < 1}, where a runs through all 
vectors in { — 1, 1}". For a = {a,a, . . . ,a) E E", the sign permutation polytope V±7r{a) 
coincides with the n-cube 7„ which can be evaluated either in terms of the 2n half- 
spaces {x G E" I — « < a:^ < a, for 1 < A; < n}, or, equivalently, with respect to 
the 2" vertex set {(±q;, ±a, . . . , ±a)}. More generally, if we consider the n-tuple of real 
coefficients a = (ai, . . . , am, Oi, . . . , 0„_m) G E" consisting of A; < n distinct coefficients 
of multiplicity rrii, 1 < i < k, then the number of distinct sign permutations, i.e. 
vertices, in the set {7r(±Q;i, . . . , ±0;^, Oi, . . . , 0„_m), tt G S'„} is 

2"n' 

. ,. ^. (14) 

2"-™ mi!...mfc_i!(n-m)! ^ ^ 

Consequently, we have it that the sign permutation polytope V-i-^(a) given by 

conv{7r(±ai, . . . , ±a„, Oi, . . . , 0„_„), vr G S'„}, (15) 

has a vertex set of order 2™n!/(n2i! . . . mk^i\{n — m)\). We see again that for m = 1, the 
sign permutation polytope V±n(a), o, = (cti, . . . , am, Oi, . . . , 0„_m.) G E", m < n, describes 
the n-dimensional cross-polytope /3„, and for m = n and k = n, the sign permutation 
polytope coincides with the n-dimensional cube 7„. For a = (a, 0, ... ,0) G E"""*"^, we 
note that the n-dimensional regular simplex «„ is a permutation polytope of degree 
n + 1 defined as the convex hull of all permutations of the vector a G E""''^. 

Finally, for what follows, we derive the volume and insphere radius for the n-cross- 
polytope. 

Lemma 4.4 Let a > 0. The volume of the n-dimensional cross-polytope 
conv{7r(±a, 0, . . . , 0), vr G Sn} of edge length \/2a is (2a)"/(n)!. 

The cross-polytope /3„ is a composition of 2" convex regions with each region 
corresponding to an orthant of /3„. The convex region determined by the positive orthant 
of /3n is an n-dimensional cone of height a. In E", the volume of the cone of height h 
over an (n — 1) -dimensional base of volume B is given as Bh/n. By induction, the 
volume of the cone in the positive orthant of /3„ is a"'/n\. Consequently, the volume of 
/3„ is (2a) Vn!. 

Next we consider the inscribed sphere of the cross-polytope. 

Lemma 4.5 The insphere of the cross-polytope conv{7r{±a, 0, . . . , 0), vr G Sn}, a > 0, 
of side length \/2a has radius rm = a/y/n. 



Proof. We maximize the n-dimensional sphere Y^^=i ^i^ ~ ''"^ subject to the constraint 
Xir=i -^j ~ ^- '^^^ constraint Yll=i Xi = a represents the coordinate sum of the points 
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on the facet of the cross-polytope in the positive orthant of E". Let L{xi, X2, ■ ■ ■ , Xn) = 
Sr=i ^«^ ~ -^ (X]r=i 3;j — «) . By the method of Lagrange muhiphers, we have it that 
VL(xi, X2, . . . , Xn) = {2xj — X,i = 1, . . . , n} with J2^=i ^i = '^- Solving the system of 
equations yields Xi = a/n for i = 1, . . . ,n, and the result follows. D 

5. Quantum sign permutation polytopes 

A fundamental topic of quantum information theory has been the characterization of 
entanglement amongst states of a quantum system. A natural question concerns the 
partitioning of the set of all states of some finite dimensional composite quantum system 
according to entanglement type. It has been shown that all bipartite entangled pure 
states are asymptotically equivalent, up to local operations and classical communication, 
to the Einstein-Podolsky-Rosen state (Bennett et al (1996)). For systems of three or 
more parties, there are several equivalence classes of different entanglement types (Vidal 
(2000)). We note the particular instance of three qubit systems in which Diir et al 
(2000) has shown that there are two equivalence classes possessing genuine tripartite 
entanglement. 

In this article, we focus on the construction of sign permutations polytopes 
within the space of quantum states, paying attention to three-qubit systems. We 
demonstrate the construction of quantum sign permutation polytopes for both pure and 
mixed quantum states, and we utilize the weak majorization feature admitted by sign 
permutation polytopes to classify a certain subset of quantum states. The construction 
of sign permutation polytopes from section |4] can be translated to the quantum world 
as follows. 

Let c/ G N and let "H denote the c?- dimensional complex Hilbert space. Let Aid 
denote the set of states that coincide with the set oi d x d complex Hermitian matrices 
with non-negative eigenvalues and unit trace which act on the Hilbert space Ti; 

Md = {g\g = g^;g>0;TTg = l}. (16) 

Aid is a (t/^ — l)-real-dimensional convex set and the vertices of this set form a 2(c/— 1)- 
dimensional subspace of the (c?^ — 2)-dimensional boundary dAid- The set of density 
matrices of A^^^ consists of pure and mixed states and the pure states correspond to 
the extreme points of this set. For g G Aid-, f? is a pure state if and only ii g^ = g 
with g written as a rank one projector onto l-i. Otherwise, ^ is a mixed state, and we 
write ^ as a convex combination of pure states. Noting that any general n-dimensional 
affine space A is isomorphic the n-dimensional afiine space E" (Br0ndsted (1983)), we 
put n = d^ — 1, and introduce a correspondence between a state ga & Aid and a point 
a G E*^'-^ by 
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Qa 



0:2 — OL^L 



Mr. 



a2 + a^L . . . \ 



-)■ 



04 



v 



H-> 



/ 



^2 
^4 



(17) 



V / 

To ensure that n- dimensional polytopes have non-zero volume in E", we require that 
the dimension of the corresponding affine hull be n + 1. This fact may be demonstrated 
by noting that affine independence of an {n + l)-family of points in E" is established if 
and only if there is linear independence amongst an augmentation of the (n + l)-family 
of points in E""*"^. 

Lemma 5.1 (Br0ndsted (1983)) Let {ai] G E", i = 0, . . . , n, and let T{a) : E" -^ E"+i 
be defined by T{a) = (a, 1). The set {ai}, i = 0, . . . ,n, is affinely independent in E" if 
and only if {{ai, 1)}, i = 0, . . . ,n, is linearly independent in K"'^^ . 

Proof. Suppose that {a^}, i = 0, . . . ,n, is affinely independent in E", i.e., ^"=0 -^^^j 
vanishes in E" and Y17=o ^i vanishes in M only when Xi = 0, i = 0, . . . ,n. In particular, 
we have (I]"=o -^^^^ Sr=o •^«) = Z]"=o '^i('^i' 1) = (Oe",Oir) = OEn+i only when A^ = 0, 
i = 0, . . . ,n. Consequently, {{ai, 1)}, i = 0, . . . ,n, is linearly independent in E"+^. 
Conversely, suppose that {(a,, 1)}, i = 0, . . . ,n, 
Yh=i f^Mi^ 1) = Oe"+i only when /i^ = 0, z = 0, . 
Y17=o l^i ~ Ok occurs only when /ij = 0, i = 0, . . . 
independent in E", and the result follows. 



is linearly independent in E"+^, i.e., 
. . ,n. Therefore, Yll=o f^i'^i — Oe" and 
n. Thus, {ai}, i = 0, . . . ,n, is affinely 

D 



5.1. Pure state polytopes 

Definition 5.2 A pure state polytope is a polytope whose vertex description is given 
by a set of pure states, i.e., conv{gi} where Qi G M.d such that Tigf = l,Qi = g], Qi > 
and TiQi = 1. ^4 pure state sign permutation polytope is a polytope whose vertex set is 
constructed by applying the set of sign permutations on a pure state and admitting only 
those sign permutations that yield a valid pure quantum state. 

Let f) G A^2 be a pure qubit state. 



Q 




I -P. 



with Px^ + Py'^ + Pz^ = 1 and P^, Py, Pz G M. The parametrization oi g E Ai2 is given by 
{P^,Py,Pz) G M^. Applying the set of sign permutations to the qubit parametrization 
(1/v^, 0, 1/v^) and taking only those sign permutations that return valid pure quan- 
tum states yields the cuboctahedron of Figure |2} 
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Figure 2. Qubit Cuboctahcdron: sign permutation polytope with basis vector 
(l/\/2, 0, l/\/2). The cuboctahedron consists of 8 triangular faces, six square faces, 24 
edges and 12 vertices. 



For the construction of a quantum sign permutation polytope from a pure state 
\ip) G "H, we enumerate the set of sign permutations on the Euchdean representative of 
density matrix \ip) {ip\ G Aid and consider only those sign permutations ivTpuro for which 
iTTpured'?/') (V'l) IS again a pure state. For three-qubit systems, the general description 
of a pure state \ip) is given by 

\^) = Ao |000) + Aie'^ |100) + A2 |101) + A3 |110) + A4 |111) (19) 

where Ai > 0, < ^ < vr and ^ . Af = 1 (Acm et al (2001)). Diir et al (2000) have shown 
that non-trivial tripartite entanglement of three qubit systems can be described in terms 
of two equivalence classes. In particular, a state in the system of three qubits possessing 
genuine tripartite entanglement may be a Greenberger-Horne-Zeilinger (GHZ) state for 
which a typical representative state is given by 

|GHZ) = 1/V2(|000) + |111)). (20) 

On the other hand, a three-qubit state may be in the W class of tripartite entangled 



states. A quantum state is said to be W-type entangled if A4 = 6^ = in (19). A typical 
representative state for those states possessing W-type entanglement is given by 

|Vr) = 1/^3(1100) + |010) + |001)). (21) 

As an example, we construct a quantum polytope of high dimensionality for 
some fixed affine space that possesses W-type entanglement. For such a polytope, we 
considered a random three-qubit pure W-state \W) = 0.7586 10) +(0.809 - 0.588;.) |2) 
+(0.809 + 0.588;) |5) +0.242 |7). Taking the set of all possible sign permutations on 
the parameterized vector, we construct the polytope vertex set by only including those 
sign permutations that correspond that valid W-type entangled pure quantum states. 
W-type entanglement is ensured by showing a vanishing 3-tangle on the set of valid 
quantum states. The 3-tangle, r^, is the measure most often utilized to distinguish 
between the GHZ- and W- class of three-qubit states (Coffman et al (2000)). For 
pure states, the 3-tangle is equivalent to the absolute value of a quantity known as 
Cayley's hyperdeterminant (Cayley (1845)), and is known to vanish for those states 
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possessing W-type entanglement. In our example, there are 2^8!/4! = 26880 possible 
sign permutations for a general three-qubit W-type entangled state of which a total of 
5376 pure points are returned as possessing W-type entanglement. 

5.2. Mixed state polytopes 

Definition 5.3 A mixed state polytope is a polytope whose vertex description is given 
by a set of mixed states, i.e., conv{gi} where Qi G Aid such that Tigf <1, Qi = q\, f?j > 
and TiQi = 1. A mixed state sign permutation polytope is a polytope whose vertex set 
is constructed by applying the set of sign permutations on a mixed state and admitting 
only those sign permutations that yield a valid quantum state. 

An open problem of multipartite entanglement is the identification of the 
entanglement type for a given quantum state. Appealing to the theory of convex sets 
may help resolve the issue of multipartite entanglement identification whereby a state 
of unknown entanglement type is described as a convex combination of states of a 
fixed entanglement type. As noise has always to be considered for an experiment, we 
are interested in describing a volume (an e-ball) around the given state which is still 
of the same entanglement type. However, achieving a lower bound on the allowed e- 
environment requires knowledge about the facet set of a polytope of a fixed entanglement 
type. For n-dimensional convex polytopes consisting of a vertex set of order m, the set 
of hyperplanes can be of the order m^^/^J (Matousek (2002)). In this instance, an e-ball 
argument to identify multipartite entanglement is rendered inefficient. 

Let us now suppose that there exists a convex decomposition of a given state 
g G Aid in terms of states of a fixed entanglement type (Kampermann et al (2010)). 
We consider the construction of quantum sign permutation polytopes whose vertex set 
is now contained in the convex hull of the given convex decomposition. The volume 
of the constructed sign permutation polytope allows conclusions about the robustness 
of state identification. For efficiency in high dimensions, we consider the restriction 
of a general sign permutation polytope to that of the n-cross-polytope. The quantum 
polytope analysis proceeds in two parts. The initial part of the process is concerned 
with the construction of the regular n-cross-polytope. In a second algorithm we address 
the robustness of the identification process. 



Algorithm 5^ constructs a sign permutation polytope of edge length \/2a, a G M, 
centered on a given state g G Aid- The resulting polytope describes a convex hull of 
states of a certain entanglement type by ensuring that the vertices of the constructed 
quantum polytope be written as a convex composition of points from the initial convex 
set of states possessing a fixed entanglement type. We assume that the convex set input 



of algorithm 5.4 has order m > d — 1. This is because we wish construct polytopes 



with non-zero volume centered on a mixed state in Aid- 
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Figure 3. Constructing a quantum (d^ — l)-cross-polytope: V±^(o), tt g 5'rf2_i with 
basis vector a = (a, 0, . . . , 0). 1. illustrates the given state gi as a convex combination of 
states of a fixed entanglement type. 2. 3. and 4. illustrate the construction of V±^(a) 
centered on q while 5. shows a (d^ — l)-quantum-cross-polytope centered on the given 
state Q € Md- 



Algorithm 5.4 Quantum polytope construction 
Input: 

(i) Given state g E Add 

(a) Convex set of states of a fixed entanglement type, Qi G M.d for i = 1, . . . ,m and 
m > d'^ — 1, such that q = Yl^i ^iQii^i ^ '^'^^ 'YlT=i -^j = 1- (Kampermann et al 
(2010)) 

Output: Quantum (rf^ — l)-cross-polytope of a fixed entanglement type with edge 
length \/2a centered at g, volume (2a)'^ '^ /{^P — 1)! and an inscribed sphere of radius 



^in = olJ\J (P — 1, see lemma 4-4 
Procedure: 



and lemma 



4.5. 



(i) Let Q G W"" "" denote the Euclidean vector associated with g E Aid- Let qi G 
E"^ ~^, 2 = 1, . . . , r/i, he analogously defined. 

(a) Let Tg(a) : a — g, a G E,'^ ^^ , be an affine transformation on E,'^ ^^. Let 



conv{Tg{g. 



Q\l^ljl 



, m} be the translate of conv{gi, i = 1, . . . , m} under T^ 



(Hi) For suitable a G M, we ensure that the set of sign permutations of the vector 
(a, 0, . . . , 0) G E^^^ is contained in the convex hull ofTg(g.i), i = I, . . . ,m. 

(iv) Maximise a G M such that the set of 2{d'^ — 1) points {7r(±a,0, . . . ,0),7r G Sd^-i} 
are still contained in the convex hull ofTg{gi), i = 1, . . . ,m. 

(v) Let T-'^^a) be an affine transformation on E'^ ^^ given by T^^(a) : a + g, a gK'^ ^^ . 
Then 



conv{Tri({7r(±a, 0, . . . , 0), tt G Sd^-i})} 



(22) 
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is the {(P — 1) -dimensional cross-polytope with edge length \/2a and volume 
{2aY ^^/{cP — 1)! centered at g and contained in the convex hull conv{gi,i = 
l,...,m}. 



Step (i) of algorithm 5.4 establishes a correspondence between the set of density 
operators Aid and vectors in Euclidean space E'^ ~^, see figure [3J (1. and 2.). The 
quantum states g and gi, for i = l,...,m, are mapped to their corresponding 
representatives in E'^ ~^. Step (ii) describes the translation transformation Tg that 
acts on the set of classical vectors g and gi, for i = 1, . . . ,m. Under Tg, the vector 
^ G E'^ ~^ is mapped to the o G E'' ~^. Correspondingly, the set gi, for z = 1, . . . ,m, 
is mapped to a surrounding neighbourhood of the origin. Step (iii) then constructs 
a sign permutation polytope with edge length y/2a centered at origin for the basis 
vector a = (a,0,...,0) G E!^~\ see figure p^ (3.). Step (iv) utilizes a divide and 
conquer algorithm (O' Rourke (1998)) to ensure that the value a G M in the basis 
vector induces the polytope of maximum volume. By lemma 4.5[ the constructed sign 
permutation polytope has volume (2a)'^ ~^ /{d'^ ~ !)'• Since this polytope is defined by 
taking the sign permutations of the vector (a, 0, . . . , 0) G E'f ~^, the construction yields 
the regular (d^ — l)-dimensional cross-polytope. Step (v) implements the translation 
transformation T^^ on the set of vertices of the constructed cross-polytope. Since a 
translate of a convex polytope is again a convex polytope, a consequence of result 
X + convM = conv(a; + M) (Br0ndsted (1983)), for x G E" and M C E", it follows that 
conv{T^^({7r(±a, 0, . . . , 0), vr G Scp-i})} is a regular cross-polytope centered at g with 
edge length ^/2a and volume {2aY '^/{(^ — l)', see figure [3J (4.). Since any (d^ — 1)- 
dimensional affine space is isomorphic to the Euclidean space E'^ ~^, the constructed 
polytope yields a desired quantum polytope in A^,^ with respect to the Hilbert-Schmidt 
norm, see figure [3] (5.). D 

The next stage in the process addresses the question of robustness of multipartite 



entanglement identification. Algorithm 5.6 makes use of the weak majorization criterion 



to determine the points g' G M.d in the convex hull of a constructed sign permutation 
polytope. The set of points g' G M.d can be represented as the state g G M.d plus 
noise in the quantum system. We make use of the following result from Zyczkowski and 
Sommers (2003). 

Theorem 5.5 (Zyczkowski and Sommers (2003)) The volume of the [d"^ — 1)- 
dimensional set of states Aid with respect to the Hilbert-Schmidt measure is given by 

Vo1hs(-M.) = v^ 2(,_,v, r(rf^) (23) 

where r(-) is the Euler gamma function. 

Algorithm 5.6 On the robustness of state identification 
Input: 

(i) Cross polytope conf { 7r(±a, 0, ..., 0), vr G Sd^-i} of edge length \f2a and volume 
(2a)'^ ~^ /d^ — 1! associated with a quantum polytope of a fixed entanglement type. 
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Figure 4. Translate of an arbitrary state g' , Tg{g'), contained in the convex hull of 
the {(P — l)-cross-polytope V±7r((j), n € S^^-i with basis vector a = (a, 0, ... , 0). 



(a) Given state g E Aid- 
(Hi) An arbitrary state g' G Aid- 

Output: The fraction of the volume of Aid describing the set of all states in the 
constructed sign permutation polytope possessing an entanglement type comparable to 
the state p G Aid- 
Procedure: 

(i) Let g' G E,^ ~^ denote the Euclidean vector associated with an arbitrary state 

g' e Md- 

(ii) Let Tg{g') be the image of g' under the mapping Tg{a) : a — g, a G E"' ^^. 

(Hi) Consider the positive cone E^^^ := {(ai, . . . ,ad2_i) G E'^ ~^\cii > 0}. For suitable 
sign changes, let T^g{g') be the image ofTg(g') in the positive cone E^~^. 

(iv) Implement corollary 4-i' Letting V-i-7r(a) denote the cross-polytope conv{'K{±a, 0, . . . , 0), 
TV G Sd2-i} of edge length \/2a. Then 

T+-,{g') G V±,(„) ^ T+-,{g') ^^ a. (24) 

(v) If T^g{g') -<w a then T+g{g') is contained in sign permutation polytope V±,r(a)- 
Consequently, T^g{g') represents a state in the 

-J (25) 

,/d7T<ii<i-^)/^r {!)... r{d) 

fraction of all states contained in sign permutation polytope V±n{a) which possesses 
an entanglement type comparable to state p E Aid- 



Step (i) of algorithm 5.6 considers a given arbitrary state g' G Aid to describe the set 



of states differing from g G Aid- By analyzing the classical representation of g' G Aid, 
a characterization of quantum states possessing an entanglement type comparable to 
V±7r(a) can be achieved. Step (i) begins by evaluating the classical representation 
of ^'; g' G E*^ ^^. Step (ii) implements the translation transformation Tg(a) on 
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g' E E"' ~^ thus translating g' E E,'^ ~^ to a neighbourhood of the cross-polytope 
V±7r(a) = conv{7r(±a;, 0, . . . , 0), vr E Sd^^i}. Step (iii) prepares the point g' E E'^ ~^ 
for the weak majorization criterion of corollary 4.1 by placing g' E E*^ ~^ in the positive 
cone E_,_ ~^ for suitable sign changes. Step (iv) then implements the weak majorization 
criterion for point inclusion in a sign permutation polytope. If T+^(^') -<^ a then 
noting Tg{g') -<w T+g{g'), it follows by the transitive property of weak majorization 
that Tg{g') -<^ a, see figure kl Consequently, Tg{g') E V±,r(a)- Finally, step (v) returns 
the fraction of states of M.^ contained in the quantum cross-polytope V±7r(a) centered 
on the state g E Aid- This fraction is given as ratio of the volume of the quantum 
cross-polytope V±,r(a) to the volume of the set of all density matrices (Zyczkowski and 
Sommers (2003)). D 

Algorithm |5.6 introduced the fraction of all states of M.d contained in a 



quantum sign permutation polytope possessing an entanglement type equivalent to the 
entanglement type of the state centered at the given quantum polytope. This fraction 
was established in terms of the ratio of Hilbert-Schmidt volumes. However, a Euclidean- 
e-ball argument describing the distance between a given state and a set of states of 



comparable entanglement type may also be considered. By lemma 4.5 , the insphere of 
the constructed quantum polytope conv{7r(±a, 0, . . . , 0), vr E Sd'i-i} given in algorithm 
is of radius a/ y/d? — 1. Letting e = a/\/d? — 1, the set of states g' E Aid contained 



5.4 



in e-ball possessing an entanglement type comparable to that of the state g E Aid is 
given by the Hilbert-Schmidt distance 



\\g' -Q\\us<(^/Vd^-l. (26) 

Now we note that the volume of an n-dimensional sphere of radius r is V^ = S'„_ir"/n 
with Sn-i = 27r"/^/r(|) denoting the hyper-surface area of an n-sphere of unit radius, 
and r(-) denoting the Gamma function (Bengtsson and Zyczkowski (2006)). For radius 
i"in = a/\/(P — 1, the {(P — l)-dimensional insphere of the cross-polytope of edge length 
\/2a has volume 

On the other hand, the volume of the (d^ — l)-dimensional cross-polytope of edge 
length \/2a is {2aY ~^ /{d^ — 1)!- Comparing these volumes yields an estimate that the 
(d^ — l)-dimensional insphere of radius otj ^d? — 1 is approximately {jx j^^'^ -i)/2 ^j^at 
of its respective cross-polytope. Consequently, while the characterization of arbitrary 
states with respect to polytope construction herein which possess an entanglement type 
comparable to the given state g E Aid is readily obtainable by both weak majorization 
and Euclidean norm calculations, arguing robustness of identification in terms of volume 
ratios between quantum sign permutation polytopes and the space of density matrices 
may be of some merit. 

Finally, the lower bound on the fraction of states introduced in algorithm 5.6| may 



be improved by considering the construction of the sign permutation polytope V±.jr{a) 
with a = («!, . . . , am, Oi, . . . , 0rf2_i„„), m < d"^ — 1, denoting a multiset of k distinct 
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elements of multiplicity rrii {1 < i < k) . This approach is more costly and requires that 
2'"n!/(mi! . . . mk_i\{n — -m)!) vertices be evaluated as convex combinations of an initially 
given convex set of density matrices with a fixed entanglement type. Furthermore, 
the resulting polytope structure loses regularity exhibited by the n-cube and n-cross- 
polytope thereby making volume calculation difficult. 

6. Conclusion 

We discussed the construction of n-dimensional polytopes of degree n in the n- 
dimensional Euclidean space E" called sign permutation polytopes, and presented a 
necessary and sufficient condition that characterizes the points of these polytopes. We 
related the construction of sign permutation polytopes in Euclidean space to the space 
of density matrices. As an application of such quantum polytopes, we considered 
the robustness of identifying those states of a fixed entanglement type. This process 
proceeded in two stages, the first stage concerned the construction of a convex polytope 
with a vertex set of a fixed type entanglement. The second stage characterized via 
weak majorization the fraction of all quantum states contained in the convex hull of 
a constructed convex polytope which possess an entanglement type comparable to the 
polytope's center state. 
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Appendix A 



Proof of lemma |3.4[ Let us suppose otherwise. Then for x E Va and y E Vb there is 
a hyperplane Hua such that 

{x,u) < a < {y,u). (A.l) 

Since T^kKuVk = Y.k<n(^k, we may take the hyperplane Hi^^j^kKnVk^ 1" = (1, • • • , 1) e 
E", to be the representative for the set of hyperplanes H^a that place V^ in one of the 
closed halfspaces of E", see Appendix B. In particular, Vb is in the halfspace H^ y 
that satisfies J2f,<^yk < ivAn)- Now, (x, 1„) = Y.k<n^k- Consider the sum Y.k<n^k- 
For X E Va there are non-negative scalars Xi,i = l,...,m with m < n\ such that 
Ylii<m -^j = 1 ^^d X = X]i<m '^i^i- Therefore, we have it that 

^ a;fc = ^ ^ XiUik = ^ ^i^ Oiik = ^ \i^ ttk, (A. 2) 

k<n k<n i<.m i<m k<n i<m k<n 

since Y.k<n'^ik = Y.k<n'^k for z = 1, . . . , m. Furthermore, since Y.k<n^k = Y.k<nPk, 



we have aj = ^ i<nPi — ^ i<nCii for j G {1, . . . , n}. It follows that equation (A. 2) may 
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(1.2,3) 



(3,1, 




(3,2,1) 

Figure Bl. Permutahedron: Va, a — (1,2,3) 



be written as 

^Aj ^fc<„afc + ^ A - ^i<„az I , (A. 3) 

i<m \ ^-^ l<n ^-^ / 

for j G {1, . . . , d}. Simplifying, we write equation (|A.3|) as X]i<m -^* "^Kn l^i- Since there 
are numbers /Xj > satisfying Yl,i<m' ^^i = 1, we then have it that Z]i<m -^^ Z]«<n /^« = 
Ei<m'/^«E«<nA- Moreover, since Efc<„/3fc = Efc<nAfc' ^^^ ^ = l,---,"^', it follows 
that Y.i<m' /"* E/<„ A may be given by 

^ /ii ^ A = 5^ 5^ /i,A = 5^ 5^ /^*A; = 5^ yi- (A.4) 

i<m' l<n l<n i<m' l<n i<m' l<n 

Therefore, (x, 1„) = (y, l.„), and this contradicts the initial assumption that the convex 
hulls Va and Vf, can be properly separated. This completes the proof. D 

Appendix B 

We construct the hyperplane Hu,a with the property Hu,a H Va 7^ 0. 

Example 1: Constructing hyperplanes for the permutahedron Va with a = (1,2,3), see 
figure |B1[ 

Taking the points (2, 1, 3) and (2, 3, 1), we have it that the parametric equation of 
hyperplane (line) containing these points is given by Hi : (0, 2, — 2)s + (2, 3, 1), for s G E. 
For any z E Hi, we may write z componentwise as zi = 2, Z2 = 2s + 3, z^ = —2s + 1, 
s G E. Now since a hyperplane may be written in the form 

H^,a = {zeE''\{z,u) = a} (B.l) 

for u G E'^\{0}, and a G M.'^, we may rewrite Hi in terms of the Hu,a form. We put 
a = Y^ak = Q- Then 2ui + {2s+3)u2+{—2s+l)us = r. In particular, we have it that if„ g 
is constrained by 2ui +4^2 = 6 with U2 = M3. We determine the remaining hyperplanes 
in a similar fashion. Proceeding anticlockwise, H2 = (1,-1, 0)s + (3,2,1), for s G E. 
Rewriting H2 in terms of an H^q expression, the candidate vector u G E^ is constrained 
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by the equation 5ui + u^ = 6 with Ui = U2. The hyperplane H3 = (0, —1, l)s + (3, 1, 2), 
for s G E, rewriting H^, the candidate vector -u G E^ is constrained by the equation 
3-Ui + 3u2 = 6 with U2 = u^. The hyperplane H4^ = (— 2,2,0)s + (1,3,2), for s G E, 
rewriting H4, the candidate vector -u G E^ is constrained by the equation 4ui + 2^3 = 6 
with Ml = M2- The hyperplane H^ = (0, 1, — l)s + (1,3,2), for s G E, rewriting H^, the 
candidate vector m G E^ is constrained by the equation ui + 5u2 = 6 with U2 = U3. 
Finally, the hyperplane Hq = (—1, l,0)s + (1,2,3), for s G E and rewriting Hq in the 
form Hufi, the candidate vector -u G E^ is constrained by the equation 3-ui + Su^ = 6 
with Ml = U2. Note that while each evaluated hyperplane has infinitely many solutions, 
a candidate that satisfies each of the particular constraints is the vector with compo- 
nents all equal to 1. Of course, since the direction vectors in each of the parametric 
expressions are permutations of each other then a candidate vector -u G E^ which sat- 
isfies all of the above constraints is precisely that vector which, under permutation of 
indices corresponding to the direction vectors, generates a solution vector u' G E^ for 
all hyperplane -ff^'.e- The only such candidate is the vector m = (1, 1, 1). 

Example 2: Constructing hyperplanes for Va, a = (1, 2, . . . , n). 

Continuing in the manner of above, we seek a candidate hyperplane for the 
set of all hyperplanes bounding Va entirely in a closed half space. We claim the 
hyperplane given by H^ n(n+i) , 1„ = (1,...,1) G E", is such a candidate. Taking 

in, 2 

the vector a = {1,2, ... ,n) G E", the associated permutahedron Va is defined as 
the convex hull of all vertices obtained by permuting the coordinates of the vector 
a. Since the set of vertices connected by and edge to the vector a are exactly those 
permutations of vector a that differ by an adjacent transposition, we have it that the 
set of vertices connected by an edge to vertex a is the set {(1,2, . . . ,n — 2,n,n — 1), 
(1,2, ... ,n — l,n — 2,n), . . . , (2, 1, . . . ,n)}. Taking the vertex a G E" together with 
any n — 2 vertices of the connecting set, of which there are n — 1 possible choices, 
determine a (n — 2)-dimensional hyperplane. Fix the hyperplane Hi to contain the 
set of vertices connected by an edge to the vector a G E" that exclude the last vertex 
(2, 1,3, ... ,n). The parametric equation for the hyperplane Hi is written in terms of 
the directions vectors {(0, . . . , 0, 1, -1), (0, . . . , 0, 1, -1, 0), . . . , (0, 1, -1, ... , 0)}. Each 
direction vector is represented as the difference between the vertex a and an element 
of the set {(1, 2, . . . , n — 2, n, n — 1), (1, 2, . . . , n — 1, n — 2, n), . . . , (1, 3, 2, . . . , n)}. 
Rewriting Hi in the Hu,a form with a = ^^^^ — -, we have it that the candidate 

vector -u G E" is constrained by the equation ui + i ^^''^ ' — 1 j ^2 = 2 with 
U2 = U3 = . . . = Un- Next, the hyperplane H2 determined by the vertex a G E" 
and its set of connecting vertices that excludes the vertex (1, 3, 2,4, . . . ,n) has its 
parametric expression written as linear combination of the set of direction vectors 
{(0,..., 0,1,-1),(0,...,0, 1,-1,0),..., (0,0,1, -1,0..., 0), (1,-1,0,...,0)}. Writing 
the parametric equation of the hyperplane H2 in terms of the H n(n+i) form, we have the 

candidate vector -u G E" constrained by -Ui + ( ^^^^ — - — 1 j -^3 = ^-^ — - with Ui = U2 and 



Quantum sign permutation polytopes 19 

Us = U4 = . . . = Un- In a similar manner, the remaining n — 3 hyperplanes containing the 
vertex a can be formed. Solving for the vector -u G E" that determines the hyperplane 
H n(n+i) associated with Hi and H2 yields the vector 1„ = (1, . . . , 1) G E". It follows 
that H n(n+i) may be taken as the representative for the set of hyperplanes containing 

in, 2 

vertex a G E". Furthermore, since all vertices a' G E" of the permutahedron Va C E'^ 
are permutations of the vertex a, for Tia' G Sn, then the set of hyperplanes containing 
the point a' G E" can be rewritten in terms of a hyperplane H n(n+i) for some candidate 

U, 2 

vector u G E"^. For such a candidate u G E" and permutation vr^/ G Sn, H^-i,, „(n+i 



1) 



is a representative hyperplane for the set of hyperplane containing the vertex a G E". 
Consequently, H^-i,, n(n+i) = H^ n(n+i) . In particular, Ti^^iu) = 1. Since the vector 
with all components equal to 1 is invariant under permutations of the indices, it follows 

that H n(n+i) may be taken as a representative for the set of all hyperplanes bounding 

^' 2 

Va entirely in a closed halfspace with the property that H^^a H Va 7^ 0. 
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